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Abstract. In the present work, we explore efficient ways to transmit binary 
information in discrete, semi-infinite chains of coupled oscillators using the 
process of nonlinear supratransmission. A previous work showed that such 
transmission is possible and, indeed, reliable under the idealistic condition 
when weak or no damping is present. In this paper, we study a more real- 
istic case and propose the design of mechanical devices in order to avoid the 
loss of information, consisting on the linear concatenation of several such me- 
chanical systems. Our results demonstrate that the loss of information can be 
minimized or avoided using such physical structures. 



1. Introduction 

The recently discovered process of energy transmission in the forbidden band 
gap of nonlinear chains pQ has opened a wide variety of new physical applications. 
This phenomenon consists in the sudden increase in the amplitude of wave signals 
that propagate in a semi-infinite, discrete, nonlinear chain driven at its end by a 
harmonic disturbance irradiating at a frequency in the forbidden band gap, and 
it has been studied in sine-Gordon and Klein-Gordon systems [2], double sine- 
Gordon systems [3], Fermi-Pasta- Ulam systems [4], Bragg media in the nonlinear 
Kerr regime [5], and even in continuous, nonlinear, bounded media described by 
undamped sine- Gordon equations driven at one end [6]. 

From a mathematical point of view, the theory that describes the mechanism 
sustaining the phenomenon of nonlinear supratransmission has not been entirely 
unveiled in general. In fact, most of the analytical results on undamped, discrete 
systems that we posses nowadays rely on the use of approximation techniques such 
as the rotating wave approximation [3], or the study of the continuous limiting 
case [2]- Meanwhile, the study of the damped case has relied on the development 
of new perturbation methods for the weakly damped scenario [7J, and the use of 
conventional or sofisticated computational techniques [8] . 

From an experimental perspective, the phenomenon of supratransmission in non- 
linear media described by sine-Gordon equations was first observed in mechanical 
systems of coupled pcndula |3| , and applications to the design of digital amplifiers 
of ultra weak signals [§] and light detectors sensitive to very weak excitations [TP] 
have been realized recently. Further applications to optical waveguide arrays using 
the discrete nonlinear Schrodinger equation [11] , the realization of light filters [12] , 
and the propagation of binary signals in undamped or weakly damped mechanical 
chains of oscillators [13] have been recently proposed. 
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In this article, we tackle the problem of propagating binary signals in damped, 
discrete sine-Gordon systems through amplitude modulation (a rather interesting 
physical problem not studied hitherto) , with future studies of signal transmission in 
Josephson-j unction arrays in mind. We make use of analytical approximations and 
perturbation methods derived from the continuous limiting situation (equivalently, 
the strongly coupled case) in order to propose the design of a propagation device for 
binary signals in damped mechanical systems. Our numerical study will be based 
on a novel computational technique with energy- invariance properties [5], and the 
proposed designs will be based on the serial concatenation of such systems. 

Section [2] of this work is devoted to introduce the model under study, analyt- 
ical results on the exact solution of the classical sine-Gordon equation, and the 
perturbation technique used; the presentation of the computational method to be 
employed, the finite-difference scheme to estimate the continuous energy and the 
numerical prediction of the supratransmission occurrence will be the topic of study 
in Section [3l Section [4] presents the mathematical procedure to generate binary 
signals, and analyzes characteristics of the propagating signals in the discrete sys- 
tem with respect to the value of the damping coefficients and the normalized bias 
current. In Section [5j we propose an alternative for the efficient propagation of 
binary signals in discrete sine-Gordon systems based on the linear concatenation 
of semi-infinite chain either bound mechanically of energetically. The results of the 
simulations for both models are presented in this section, and we summarize our 
conclusions and propose new directions of future research in a final stage. 



2. Analysis 



2.1. Mathematical model. Throughout this paper, we assume that a, f3, 7 and c 
are nonnegative real numbers. Likewise, we consider a system {u n )'^L l of particles 
satisfying the mixed- value problem studied in [8], namely, 



(1) 



dt 2 
subject to : 



'dt 



A*U„ + P 



du n 

~~dt 

= 0. 



«n(0) 
^(0)=0, 
uo(t) = V>(*)> 



V'{u n ) = 
n e Z+, 
n e Z+, 
t > 0, 



0, 



where a and P evidently function as the internal and external damping coefficients, 
respectively. Here, A^u„ is used to denote the spatial second-difference u n+ \ — 
2u n + w„_i for every n G Z + , the boundary-driving function is given by ip(t) = 
A(t) cos(fii) for every t G (0, +00), and V(u n ) — 1 — cos(u„) — 7w„ where, due to 
its analogy with the model describing long Josephson junctions, the parameter 7 
will be denominated the normalized bias current. Notice that the Hamiltonian of 
the n-th site is 

(2) H n = - [ill + c 2 {u n+1 - u n ) 2 ] + V(u n ), 



for any differentiable function V. After including the potential energy from the 
coupling between the first two oscillators, the total energy of the system becomes 



(3) 



E 



71=1 
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Figure 1 . Bifurcation diagram of occurrence of critical amplitude 
versus driving frequency for problem ([T]) with c = 4 (solid line). 
The limiting solution to dU is depicted as a sequence of crosses. 



The next result provides us with a convenient way to compute the instantaneous 
rate of change of the energy for solutions of (JXJ) with square- summable derivatives in 
time, that is, solutions (u n (t))^ =0 for which ^ u\{t) is convergent for every t > 0. 

Lemma 1 (Maci'as-Di'az and Puri [8]). Let (u n (t))^ =0 be a solution of system 
such that (u n (t))^L is square- summable at any fixed time t. The instantaneous 
rate of change of the total energy in the system is given by 



'{uq - ui)u - /3^2(u n ) 2 - a 



^(iin - u n -i) 2 + {ui - u Q )u 

,n=l 



□ 



Corollary 2 (Geniet and Leon [3]). Let (u n (£))5£L be a solution of undamped 
system dXJ) with 7 = 0, such that (u n (t))^L is square- summable at any time t. 
Then 

(5) E(t) = -c 2 f u (s)( Ul (s) - u (s)) ds 

Jo 

□ 

We wish to mention that qualitative results similar to those presented here were 
obtained using the expressions of the individual energies provided by |14) . These 
formulas were used for validation purposes only. 
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Figure 2. Local energies H n (t) of unperturbed system (JTJ) vs. site 
n and time t, corresponding to the transmission of binary signal 

'mi', with n = 0.9. 



2.2. The continuous approximation. For theoretical purposes, we will provi- 
sionally assume that a = p = 0, and will consider a long, linear array consisting 
of N coupled oscillators described by the system of differential equations in (JTJ) , 
with a coupling coefficient 1 < c <C, N. Under these conditions, the problem under 
study can be approximated by the continuous, mixed-value problem 



(6) subject to 



d 2 u d 2 u . , 

u(x,0) = 0, 0<x<L, 
Ou 

< — (x,0)=Q, 0<x<L, 

dt ~ ~ 

u(0,t) =i/j(t), t>0, 



where L = N/c, which in turn can be simplified to the study of a sine-Gordon 
model when the normalized bias current is equal to zero. For a sufficiently long 
linear chain, a sufficiently short interval of time [0, T], and a driving function tp(t) = 
Aocos(ilt), it is reasonable to impose the boundary condition u x (L,t) = 0. 
Following |S] now, there exist solutions to problem ^ in the form 

(7) u(x,t) = 4arctan[X(z)T(t)] ! 

with X(L) — Al and normalized time scale T satisfying T(to) = 1 for some time 
t$. We assume that the solutions are periodic in time with period equal to the 
driving period, and that the solutions adapt to the driving amplitude, that is, 
Aq = 4arctan(a), which yields X(0) = a. These solutions are seen to depend on a 
free parameter A and the parameter 

(8) r = ^ 2 + A7rw 
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Figure 3. Time-dependent graph of the amplitude function used 
to generate a bit of '1' in a semi- infinite, discrete, mechanical array, 
with unitary amplification coefficient. 



The relations between a, Al and A are summarized in the following result. Here 
IK(to), sn(-,m), cn(-,m) and dn(-,m) are the complete elliptic integral of the first 
kind, the sine-amplitude and the cosine-amplitude Jacobi elliptic function of mod- 
ulus to, and the identity-amplitude Jacobi elliptic function of modulus to, respec- 
tively. 

Lemma 3 (Khomeriki and Leon [B]). 

(1) If X > then X = ALCii(k(x — L),fi) and T = cn(oj(7j — to),v), with 
equations f2K(z/) = and a = Al cn(fcL, fi) being satisfied for 

u? = A(l + T), v 2 ' 



( 9 ) / I \ 1 TA 2 

(2) If A < and TA\ < -1 then X = A L dn(fc(x - L),n) and T = sn(w(i 
t\),fjb), with equations HK.(y) = ^lu and a = Aj J dn.(kL ) fi) satisfied for 

1 

f ' 1 

1 





OJ 2 


= AT, 


V 2 


(10) 


k 2 


= ATAl, 


I? 




ti 
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TO + 
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Figure 4. Time-dependent graphs of the position and maximum 
energy amplitude of the breather generated in unperturbed system 
© by bit '1', with fi = 0.9. 



(11) 



(3) If A < and TA\ > -1 then X = A L dn _1 (fc(a; - L), fi) and T = sn(w(t 
tx), (j>), where QK.(v) = ^uj and a = Al dn _1 (fcL, /i) are satisfied for 

. n A n 



□ 



and for lo 2 , v 2 and t\ as in (2). 

Corollary 4 (Khomeriki and Leon [6]). The limiting solution of mixed-value prob 
lem ^) as L — > oo is given by 

(12) u(x,t) = 4arctan a s — — ^ LL 

\ cosh(K.Tj 

when if}(t) = Aosin(f2t). Here a 2 — A 2 (l — /i 2 ) -1 and the rest of the parameters 
are as in (3) of Theorem^ □ 

2.3. Perturbation analysis. In this section, we will assume that u(x, t) represents 
an exact solution of the unperturbed sine-Gordon equation in ([5]) with — oo < x < 
oo, and that (3 and 7 are small, nonnegative real numbers. We will consider the 
perturbed sine-Gordon model 

d 2 v d 2 v dv 

( 13 ) W~dx^ + ^ + SmV - 1 = °> 

Following }7[ , we assume the existence of a solution of the perturbed sine-Gordon 
equation of the form u(g(cj)(t))x, T(t), (j)(t)) with 



(14) 



T(t) 



9{<t>( s ))<Pi s ) ds 



and 



5(0) 



1± 
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Figure 5. Bifurcation diagram of occurrence of critical amplitude 
versus driving frequency for problem ([1]) with c = 4, for (3 = 
(solid), 0.1 (dashed), 0.2 (dashed-dotted) , 0.3 (dotted). The lim- 
iting solution to unperturbed system © is depicted as a sequence 
of crosses. 



for some suitable function </>. If we let a be a localized intrinsic mode (static 
breather) solution of the unperturbed model confined in some region |x| < 5, then 

du du 



(15) 



dt ~ dT 9 



and it follows that the Hamiltonian of the unperturbed sine-Gordon equation can 
be approximated through the expression 



(16) H*W)) 



1 f du 



2 { dT 9 



du 

dx 



- ^ +1 



cos u 



dx. 



It can be shown then that the relation between the variables 0, t and T is described 
through the system 

2 



(17) 



1 



dT 

dt 
dT 



dH* 



fduY 

Vdt) 



'7 



du 
dt 



dx, 



3. Computational approach 

3.1. Finite-difference schemes. We consider a finite system of N differential 
equations satisfying ((TJ, and a regular partition = to < t\ < ■ ■ ■ < tu = T of 
the time interval [0, T] with time step equal to At. For each k = 0, 1, . . . , M, let us 
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FIGURE 6. Local energies H n (t) of perturbed system ([T]) with ft = 
0.01 vs. lattice site n and time t, corresponding to the transmission 
of binary signal '1111', with = 0.9. 



represent the approximate solution to our problem on the n-th lattice site at time 



tk by u„. If we convey that d t u^ 



,fe+i 



,fc-l 



, that S^Un 



,fe+i 



2u» 



,fc-l 



and that <^it„ 



2u k : 



take then the discrete form 



(18) 



(Aty 



2At 



St 



the N differential equations in our problem 
V{u^)-V{u^) 



k 

x "'n 



8lu, 



2At 5t < 



k+1 



.fe-1 



0. 



where ft' includes both the effect of external damping and a simulation of an ab- 
sorbing boundary slowly increasing in magnitude on the last N — No oscillators. 
More concretely, we let un+i = un at all time, and let ft' be the sum of external 
damping and the function 



(19) 



ft"{n) = 0.5 



1 + tanh 



2n - N + N 
6 



where usually N = 50 and N > 200. 

Finite-difference scheme dTBl (which is a modified version of the one developed in 
[15] to compute radially symmetric solutions of modified Klein-Gordon equations) 
is consistent with our mixed-value problem and conditionally stable, having the 
inequality (cAi) < 1 + (a + ft' /A) At as a necessary condition for stability when 
ft' is assumed constant [8]. Moreover, if the energy of the system at the k-th time 
step is computed using the expression 



(20) 



Eh 



M 
n=l 



6 / fe+1 

yK 



fe+i 
u 
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Figure 7. Time-dependent graphs of the position and maximum 
energy amplitude of the breather generated by bit '1' in (TJJ with 
P = 0.01 and Q = 0.9. 
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where the individual energy of the n-th lattice site is computed by 
(21) K = \ ^ +X ~ + - u^)(u k n+1 - u k n ) 

V« +1 ) + V«) 
+ 2 

then the discrete rate of change of energy turns out to be a consistent approxi- 
mation of order O(At) 2 for the corresponding instantaneous rate of change, while 
the approximations of the individual and the total energies are consistent estimates 
order O(At) of their continuous counterparts. 

3.2. Supratransmission threshold. It has been proved that a discrete or con- 
tinuous medium described by damped problem (jTJ) or undamped problem ([5]), re- 
spectively, can undergo nonlinear energy transmission when the driving amplitude 
is increased above a critical threshold that depends on the driving frequency, and 
which has been numerically predicted in [6] . In the discrete case particularly, 
the numerical results obtained using the computational technique described in the 
previous section have produced bifurcation diagrams for various values of a, /3 and 
7, which are in excellent agreement with the undamped, unbiased formulation of 
the problem for large coupling coefficients. In this case the bifurcation threshold 
can be approximated via the continuous approximation provided by Corollary 21 
which in turn can be approximated through 



(22) A s — 4 arctan 
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Figure 8. Time-dependent graph of decay of a static breather in 
infinite system ((6j, for constant damping (3 = 0.01. 



For the sake of convenience, we have included a bifurcation diagram in Figure [T] 
which shows the numerical (solid line) and theoretical (sequence of crosses) critical 
amplitude at which supratransmission starts for each frequency f2 in the forbidden 
band gap. The theoretical approximation was realized using (|22j) . 

In the continuous limit, the mechanism of the bifurcation is explained by the 
fact that the medium starts to transmit energy in the form of moving breathers 
once the driving amplitude reaches its critical value [6], and the medium continues 
transmitting energy well until a vanishing driving amplitude is reached. In the 
discrete case, it has been checked numerically that this phenomenon still occurs [5] , 
and that the expression of the moving breathers is given by 



(23) 



i(i) =4 arctan 



t—vn/c 



, Vl+r 2 VT- 



cosh I 



Vl+r 2 VT- 



z/c—vt 



where r and v represent, respectively, the group amplitude and the group velocity 
of the breather. Numerical experiments have shown that r and v are variables that 
linearly depend on A s , for values of 0.87 < ft < 1 (see [5]). Moreover, we have 
verified numerically that the amplitude of the emitted breathers tend to decrease 
as the driving amplitude is decreased from A s toward zero. 



4. Numerical results 

Following |13j . in this section we consider a semi-infinite, mechanical system 
described by ([T]) . A single binary bit b will be transmitted into the medium during 
a fixed and sufficiently long period of time P equal to an integer multiple of the 
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Figure 9 . Bifurcation diagram of occurrence of critical amplitude 
versus driving frequency for problem ([lj with c = 4, for a = 
(solid), 0.1 (dashed), 0.4 (dashed-dotted) , 0.7 (dotted). The lim- 
iting solution to unperturbed system (j6]) is depicted as a sequence 
of crosses. 

driving period, by letting 

(24) A(t) = h>CA s (e-°*/ 4 - 5 - e^ - 45 ) , 

for every t € [0, P], where A s represents the amplitude threshold at which nonlinear 
supratransmission begins, and C is an amplification constant that depends of Q. It 
is worth noticing that A is identically equal to zero for a bit b equal to zero; when 
b is equal to one, the qualitative behavior of A over a period equal to several times 
the driving period is presented in Figure [3] 

Before we study the effects of damping and the normalized current in the moving 
breathers generated by the amplitude function above, it is highly convenient to 
introduce some terminology. Thus, for a breather moving away from the irradiating 
source at a constant velocity v, we assume the existence of constants <5 > and 
no G Z with the property that the breather is approximately equal to zero outside 
of the region A defined by \n — (n + vt) | < 8. In such case, we say that the position 
of the breather at time t is equal to n + vt, and that its amplitude and energy are, 
respectively, 

(25) Ab(t) = max |u n (£)| and Eb(t) = max H n (t). 

Of course, this definition of position of a breather is rather subjective. A bit more 
formally, if we know that the breather is located in a region R at time t, its position 
may be defined as the integer n\ £ R such that H ni (t) = max H n (t). 
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Figure 10. Local energies H n (t) of (JXJ) with a — 0.05 vs. lattice 
site n and time t, corresponding to the transmission of binary signal 
'1111', with ft = 0.9. 



4.1. The ideal case. Let us consider an undamped, mixed- value problem (p} for 
which the normalized bias current is equal to zero, the coupling coefficient is equal 
to 4 and the period of bit generation equals 20 driving periods. Numerically, we 
consider a finite chain system of length N = 600, a time step of 0.05, and use the 
computational scheme sketched in Section lOl A bifurcation diagram similar to the 
one found in [8] is readily obtained, and the critical amplitudes A s are used in the 
determination of the velocities for single moving breathers. Under this situation, a 
binary signal of '1111' is transmitted into chain ([1]) driven at a frequency equal to 
0.9, with an amplification constant equal to A s . 

Figure [2] shows the development of the local energy in this system with respect to 
time and lattice site. The fact that moving breathers are generated in the chain is 
clear from the graph. The fact that the first of the breathers moves at approximately 
a constant speed is shown by the upper graph in Figure |31 while its maximum local 
energy is analyzed in the bottom portion of the same figure. Clearly, this maximum 
is above a positive cutoff limit, say, equal to 5. 

In order to find a more reliable estimate of the lower cutoff limit for the maximum 
local energy carried by a breather generated using (|24|) . we have transmitted a single 
bit T' into a chain consisting of 6000 lattice sites during a period of time equal to 
2000, with all other parameters as in the previous paragraph. As a result, we obtain 
that the absolute minimum and maximum values of the maximum local energy in a 
single breather are equal to 5.6965 and 17.9835, respectively, whence it particularly 
follows that a lower cutoff limit of 5 would actually prove to be adequate. 

4.2. External damping. The next step in our investigation will be to describe 
the effects of internal and external damping, and the normalized bias current on 
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Figure 1 1 . Time-dependent graphs of the position and maximum 
energy amplitude of the breather generated by bit '1' in (pj with 
a = 0.05. 



the behavior of propagating signals in ((T|). More concretely, we wish to establish 
the qualitative effects of these parameters on the local energy amplitude and speed 
evolution of the breather solutions with respect to lattice site. To that effect, we 
consider the same mixed- value problem studied in the previous section with similar 
choices of parameters, and provide a bifurcation diagram for several values of /3 in 
Figure [5] 

The binary signal '1111' will be propagated again using the driving amplitude 
function (|24[) assuming that the medium is externally damped with a coefficient 
of 0.01. In this context, Figure [6] shows the time and lattice site evolution of the 
local energy carried by the breathers generated by the transmitted signals. These 
results evidence a drastic decay in the amplitudes of the local energies of a breather 
as it moves in the chain, a claim that is confirmed by Figure [7j in which we have 
tracked the position and the development of the maximum local energy of the 
moving breather generated by a single bit '1'. 

Two conclusions are now evident from the graphs in Figure [71 first, the breathers 
generated using our transmission process move away from the source at a constant 
speed; second, a decrease of the maximum local energy amplitude carried by the 
breather with respect to time (and, hence, with respect to lattice site) is now 
apparent. 

From Section 12.21 the solution to problem (|6]) in the continuous limit case is 
given by 

,„„n , x / Acsin(fW) \ 

26 u(x, t) = 4 arctan ^—^ , 

\l2cosh(A:r) ) 
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Figure 12. Bifurcation diagram of occurrence of critical ampli- 
tude versus driving frequency for problem ([1} with c = 4, for 7 = 
(solid), 0.01 (dashed), 0.02 (dashed-dotted) , 0.03 (dotted). 

where A = arccosh (l + (1 — J7 2 )/2c 2 ) . Using this solution and applying the pertur- 
bation method of Section [2~3l with g((f>) = 1/yl + (j> 2 , we obtain the time evolution 
of the amplitude at the origin of a static breather initially described by (|26|). with 
an external damping coefficient equal to 0.01. The results, which are presented 
in Figure |8l evidence a good agreement with the time evolution of the amplitude 
decay of the breather presented in Figure [5] 

4.3. Internal damping. Consider again problem ([1]) with the same parameters as 
in Section 14.21 and similar numerical setting, describing a medium with no external 
damping and internal damping coefficient equal to 0.05. Again, the binary signal 
'1111' will be propagated into the medium using driving amplitude function (|24|) . 
For convenience, Figure |9] provides a bifurcation diagram of energy transmission for 
several values of a. 

Under this state of matters, Figure [10] presents the dependence of the local 
energy amplitude of each lattice site with respect to time. The graphs shows a 
decrease in time of the local energy amplitude of the generated breathers. To verify 
this claim, Figure [TT] presents the time development of the position and maximum 
local energy attained by the first breather produced. This last figure also suggests 
that the breathers generated in this case move away from the driving boundary 
with approximately a constant speed, a fact that was corroborated using linear 
regression analysis. 

We must remark that, like in the case of external damping, higher values of 
internal damping produce graphs of amplitude and local energy that decrease faster 
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FIGURE 13. Local energies H n (t) of (P) with 7 = 0.1 vs. lattice 
site n and time t, corresponding to the transmission of binary signal 

'mi', with n = 0.9. 



toward zero in time, as it is expected. In any case, we have verified that the phase 
velocity of the emitted breathers remains constant. 



4.4. Normalized bias current. Finally, we carry out a similar analysis for the 
normalized bias current. 

Figure [12] shows the effect of several values of the parameter 7 in the occurrence 
of the supratransmission threshold. As before, the binary signal '1111' is generated 
by the driving boundary of a mechanical chain with parameters as in the previous 
section, and 7 = 0.1. Figure [TBI shows then the motion of the generated breathers in 
our mechanical system, while Figure Q3] provides time-dependent graphs of velocity, 
amplitude and local energy of the breather generated by the first bit. For this 
choice of 7, a slight decrease in the amplitude and the local energy of the breathers 
is shown. 



5. Application 

5.1. Concatenated mechanical chains. Let t be a positive integer, and for every 
j = 0, 1, ...£, assume that flj and Cj are positive numbers, that ctj, j3j and jj are 
nonnegative real numbers, and that Lj is a positive integer. Suppose that the 
sequence (u„)J° =1 of vectors u„ = (uj n )j—i satisfies the system of concatenated 
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Figure 14. Time-dependent graphs of the position and maximum 
energy amplitude of the breather generated by bit T' in (pj with 
7 = 0.1 and ft = 0.9. 



equations 



(27) 



dt 2 



d 
3 dt 



A 



subject to 



for j =0,1,. 



d.u 



a dllj n 
' Pj dt 
3 n(0)=0, 



dt ®=°> 
Uj (t) = 1pj(t), 



V'(u jn ) 
n e Z+, 

n e z+, 
t > 0, 



= o, 



with potential functions given by Vj(u) = 1 — cos(it) — ~1jU, and boundary-driving 
functions ipj(t) = Aj(t) cos(Vljt) for every t e (0, +oo) and j = 0, 1, . . . ,£, where 
each frequency flj is assumed to be in the forbidden band gap of the j-th discrete 
array. The number £ will be called the length of the system of concatenated arrays. 

Let Pq be equal to a multiple of the period of the harmonic driving in the zeroth 
mechanical array. As in Section 21 a single bit of information b will be transmitted 
during a period of time equal to Pq by modulating the amplitude of the driving 
oscillator in the zeroth array through 



(28) 



A (t) = 



-bCA s 
5 



0, 



-fi t/4-5 



-n t/0.45 



if b = 1. 



if b = 0. 



The reception device of the first mechanical array will be initialized then at time 
Lq/v, where v represents the phase velocity of the generated breathers in the zeroth 
array. After initialization, a time equal to a period of the driving oscillator in the 
first array will be spent in order to determine if a bit equal to '0' or T' passes by 
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Figure 15. Time-dependent local energies at sites 300 in a linear 
concatenation of mechanical arrays of length 5, as a result of the 
transmission of signal '10111001011011101001' and using a cutoff 
limit equal to 2. The decoding of the transmitted message is shown 
in each array at the top of the reception periods. 



site Lq. This determination is done by comparing the local-energy amplitude of 
site Lq with respect to a predetermined cutoff limit. Next, a driving amplitude in 
the first array defined by (|28|) will be used for a period of time equal to P\ , and the 
process continues inductively. 
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Following the conventions of SectionOTTJ the j-th system of differential equations 
in model (|27[) will be approximated numerically by a finite number of equations 
Nj Lj via the set of finite-difference schemes 

[ ' (At) 2 V J + 2At V X 3n + 2At * ° n + ~ U ' 

where /3j is the sum of the external damping coefficient in the j-th mechanical array 
and the effect of an absorbing boundary in the last few lattice sites. We mimic now 
the rest of Section [3. f I in order to define computational schemes to approximate the 
total energy and the individual energies in each mechanical array. 

5.2. Simulation. As an application, assume that the binary signal '1 Of f f OOf Of f 0- 
11101001' will be transmitted in a linearly concatenated, mechanical system (|2"T|) 
of length 5, in which a period of single-bit generation equal to 20 driving periods 
will be used. For the sake of simplicity, we will assume that otj = 0.02, j3j — 0.003, 
7j = 0.01 and Qj = 0.9 for every j = 0, 1, . . . , 5, for which numerical simulations 
yield a common critical amplitude A s = 1.778 and a phase velocity equal to 2.395. 

Computationally, we consider arrays consisting of N = 600 sites each, and as- 
sume that the link between two consecutive arrays is invariably placed on site 
L = 300. An analysis similar to that done in Section [4] reflects that the estimate for 
the minimal value of the energy carried by a breather at site L is above the value 
of A s when an amplification coefficient C — A s is used. 

In this state of matters, Figure [TBI presents the time-dependent local energy of the 
300-th site in each of the 6 mechanical arrays in our system, as a consequence of the 
transmission of signal '10111001011011101001'. For the sake of convenience, each 
graph shows approximate periods of reception and, the value of the transmitted 
bit on top of each period. Our numerical results establish that a rather reliable 
transmission can be achieved using model (|27|) under suitable parametric conditions. 

6. Conclusions and perspectives 

In this work, we have presented a simple, mathematical model to efficiently 
transmit binary signals in certain semi-infinite, discrete mechanical systems using 
the process of nonlinear supratransmission. The model is physically described as a 
concatenation of a finite sequence of n such arrays, in such a way that the amplitude 
of the driving oscillator at the end of the A:-th lattice for 2 < k < n, is equal to 
the local energy amplitude of some fixed oscillator in the (k — l)-st lattice, located 
relatively far away from its driving oscillator. Of course, this energy-bound series 
of arrays requires that the local energy of the linking site in the (k — l)-st lattice 
induced by the moving breather be at least equal to the critical amplitude of the 
driving disturbance in the k-th lattice, and simulations of this type of structures 
have shown that perfect transmission may be achieved under suitable choices of 
parameters. 

As a means in order to procure a better transmission, this work presents bifurca- 
tion diagrams of occurrence of nonlinear supratransmission in non-ideal situations. 
Thus, numerical predictions of this process have been provided in the presence of 
external and internal damping, and nonzero values of what we have called normal- 
ized bias currents. In the former delay in the appearance of the critical 
amplitude at which supratransmission starts is observed, as well as a decrease in 



PROPAGATION OF BINARY SIGNALS 



19 



the region where harmonic phonon quenching takes place. In the case of external 
damping, an expansion in the forbidden band gap region and a delay on the occur- 
rence of the supratransmission threshold are clearly noticed. On the other hand, 
the normalized bias current tends to rush the appearance of supratransmission. 
Likewise, our results show that the localized structure produced in our model move 
away from the boundary at constant speeds, even in the presence of external or 
internal damping, or nonzero normalized bias current. 

Before closing this paper, it is important to recall that our mathematical model 
(which is, after all, a Dirichlet boundary- value problem with zero initial conditions) 
propagates localized energy in the form of moving breathers. Naturally, question 
arises as to how different the propagation of energy will be in the case of a Neumann 
boundary-value problem, where one immediately recognizes a system of Josephson 
junctions coupled through superconducting wires. In this case, the existence of 
a nonlinear supratransmission threshold has been analytically and numerically es- 
tablished for an undamped system without normalized bias current |16) . Moreover, 
the existence of a bifurcation value under which energy transmission is drastically 
suspended — a process we have called nonlinear infratransmission in an upcoming 
paper — has been reported in the same work; however, no bifurcation analysis has 
been given for it yet, not even in the ideal case. 

In view of these facts, it is particularly interesting to describe numerically the 
occurrence of the nonlinear supratransmission and infratransmission thresholds un- 
der the presence of damping and nonzero normalized bias currents. Moreover, it is 
indispensable to determine the nature of the localized structures transmitted into 
discrete Josephson-junction arrays by the driving boundary. Preliminary results 
have shown that energy is transmitted into these systems in the form of kinks and 
anti-kinks, the properties of which are still important topics of study. 

Another interesting direction of investigation arises from the fact that certain 
type of nonlinear differential equations (amongst which is the sine-Gordon model 
[17] ) admit localized solutions that stabilize by dissipative effects [18]. These types 
of solutions have been studied in continuous models, thus the existence solutions 
in a discrete system of coupled sine-Gordon equations would also be an interesting 
question to attack. 
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